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ABSTRACT 

We report on initial results from a campaign to obtain optical imaging of a sample of Ultra Compact 
High Velocity Clouds (UCHVCs) discovered by the ALFALFA neutral hydrogen (HI) survey. UCHVCs 
are ALFALFA sources with velocities and sizes consistent with their being low-mass dwarf galaxies 
in the Local Volume, but without identified optical counterparts in existing catalogs. We are using 
the WIYN 3.5-m telescope and pODI camera to image these objects and search for the presence of 
an associated stellar population. In this paper, we present our observational strategy and method for 
searching for resolved stellar counterparts to the UCHVCs. We combine careful photometric measure¬ 
ments, a color-magnitude filter, and spatial smoothing techniques to search for stellar overdensities in 
the g- and Uband pODI images. We also run statistical tests to quantify the likelihood that whatever 
overdensities we find are real and not chance superpositions of sources. We demonstrate the method 
by applying it to two specific data sets: WIYN imaging of Leo P, a UCHVC discovered by ALFALFA 
and subsequently shown to be a low-mass star-forming dwarf galaxy in the Local Volume; and WIYN 
imaging of AGC198606, an HI cloud identified by ALFALFA that is near in position and velocity to 
the Local Group dwarf galaxy Leo T. Applying the search method to the Leo P imaging data yields an 
unambiguous detection (>99% confidence) of this galaxy’s stellar population. Applying our method 
to the AGG198606 imaging yields a possible detection (92% confidence) of an optical counterpart 
located ~ 2.5 arc minutes away from the centroid of AGG198606’s HI distribution but still within 
the HI disk. We use the optical data to estimate a distance to the stellar counterpart between 373 
and 393 kpc, with an absolute magnitude of Mi = —4.67 ± 0.09. Combining the WIYN data with 
our previous estimate of the HI mass of AGC198606 from WSRT imaging yields an Hl-to-stellar mass 
ratio of ^45—110. 
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1. INTRODUCTION 

A current issue in cosmology is the disagreement be¬ 
tween models and observations regarding the number and 
mass distribution of low-mass galaxies that exist in en¬ 
vironments like the Local Group. Simulations of the for¬ 
mation of structures like the Local Group in a ACDM 
universe predict large numbers of low-mass da rk matter 
halos around the Milky Way and M31 (e.g., Moore et 


al. ||2006[ IDiemand et al. ||2007| IGarrison-Kimmel et al. 


201 4p. Observational campaigns aimed at searching for 

low-mass galaxies in our local neighborhood, while suc- 



Willman et al. 120051 

Martin et al. 12013 

The DES 

Golla 

boration et af.|2015 ) h 

ave nevertheless not yet found 


the so-called “missing satellites” problem (Kauff mann et 
al. 1993 Klypin et al. ||1999 Moore et al. ||l999 r 


theoretical results can be mitigated when internal evolu¬ 
tionary processes (e.g., “feedback” from star formation 
and supernovae) and external processes (e.g., reioniza¬ 
tion, tidal stripping, ram pressure stripping) are taken 
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into account. The details of how and when these pro¬ 
cesses occur, and what effect they ha.ve on the galaxies, 


are complicated (e.g., Bullock et al. 2000 Guo et al. 


2010[ |Papastergis et aT 
ever, may be that the co' 


. || 2012 |. 

iold baryc 


The net effect, how- 
baryons in nascent galaxies 
become so depleted that it is impossible for the galax¬ 
ies to form stars. Alternatively, feedback may limit the 
star formation that does occur to only one generation 
of stars, rendering any stellar population that is present 
undetectable given our current observational capabilities 
and the limitations of the surveys completed so far. 

In order to both explore the “missing satellites” prob¬ 
lem and to investigate the complex interplay between 
low-mass dark matter halos, their baryonic components, 
and the gas and stellar processes that shape them, we 
have initiated a project to obtain optical observations 
of a specific class of candidate low-mass galaxies discov¬ 
ered via their neutral hydrogen emission. These objects 
were identified by the Arecibo Legacy Fast ALFA (AL¬ 
FALF A) neutral hydrogen (HI) survey (|Giovanelli et al. 
'2005|). The ALFALFA survey uses data from the Arecibo 


305-m telescope to map 7000 deg^ of the sky in the HI 21 
cm line at an angular resolution of ~ 3L5 and a spectral 
resolution of ~5 km s“*. Data acquisition for the AL- 
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FALFA survey was completed in October 2012; full pro¬ 
cessing of all the survey data cubes has been concluded, 
and catalogs have been produced for all sources at cz < 
3000 km s“^. ALFALFA has detected a sample of iso- 
lated ultra-compact high-velocity HI clo uds (UCHVCs; 
Giovanelli et al. 2010 Adams et al. 20131 with kinematic 
properties that make them likely to be located within the 
Local Volume, but not part of the Milky Way. The se¬ 
lection of these UCHVCs is described in lAdams et al. 1 
(2013); briefly, sources included in the sample must have 
kiSfll > 120 km s“^ to exclude Galactic HI sources, a 
diameter less than 30' in size, and a signal-to-noise ra¬ 
tio > 8. Sources must also be well-isolated from larger 
HI structures and previously known high-velocity clouds. 
Notably, while they are compact, the UCHVCs are gen¬ 
erally spatially resolved. The members of the UCHVC 
sample possess HI properties similar to those of Leo T 


(separations of 1°.2 and 16 km s ^ , resp ectively), earn¬ 
ing it t he nickname “Friend of Leo T”. Adams et al. 


(Irwin et al. 2007), which is the faintest known Local 


Croup dwart galaxy with evidence of recent star forma¬ 
tion {My ~ —7). At a distance of 1 Mpc, the UCHVCs 
would have HI masses of ^ 10^ — 10® Mq, HI diame¬ 
ters of ~ 2 — 3 kpc, and dynamical masses of ~ 10^ — 10® 
Mq , placing them in the range of lower mass dark matter 
halos. We are using the WIYN ObservatorjQ to image 
the UCHVCs in order to investigate whether these ob¬ 
jects are indeed nearby dwarf galaxies and, if possible, to 
study their optical emission and stellar populations. 

A dwarf galaxy located just outside the Local Group 
was recently identified using this type of approach: Leo 
P, a gas-rich, star-forming dwarf galaxy with My ~ —9.4 
and a distance of ~ 1.7 Mpc, was discovered via its 
ALFALFA HI detection and follow-up optical imagin g 
with WIYN (Giovanelli et al. |2013 Rhode et al. 12013 ). 
Several other studies have employed si milar strategies 
to loo k for nearby dwarfs. For example, Tollerud et al.| 
(2015) targeted HI sources found in the GALFA HI sur- 


vey and identified two objects in WIYN follo w-up imag¬ 
ing that ma y be Local Volume dwarf galaxies. |Bellazzim 


et al. (2015) used deep imaging with the Large Binocular 


(2015a) used the Westerbork Synthesis Radio Telescope 
(WSRT) to obtain follow-up observations of AGC198606 
and derived an HI mass of 3.5 x 10®M©, an intrinsic 
velocity dispersion of 9.3 km s“^ , a rotation velocity of 
^ 14 km s“^ , and a dynamical mass of 3.5 x lO^dMpc 
Mq. 

In this paper we present our observational strat¬ 
egy, analysis methods, and techniques for detecting 
stellar populations in the optical follow-up images of 
the ALFALFA-detected UCHVCs. We then demon¬ 
strate these detection methods by showing results from 
WIYN imaging of two fields: AGC208583 (Leo P) and 
AGC198606 (Friend of Leo T). Not surprisingly, apply¬ 
ing the method to the former held yields an unambigu¬ 
ous detection of the resolved stellar population of Leo P. 
Applying the method to the AGC198606 held yields a 
possible detection of an optical counterpart to this HI 
source, even though (as noted in Adams et al. |201^ ) no 
conspicuous stellar counterpart is present m the images. 

The paper is organized as follows: in Section]^ we de¬ 
scribe the WIYN observations, image processing, source 
detection, and photometric measurements. Section 
presents the hltering and smoothing steps that are per¬ 
formed for each set of images to search for optical coun¬ 
terparts to the ALFALFA sources. Section shows 
the results of this process for the Leo P held and the 
AGC198606 held. The statistical signihcance and the 
implications of the possible detection of a counterpart 
to AGC198606 are also discussed. The last section sum¬ 
marizes our main conclusions and describes the current 
status of our optical follow-up imaging campaign. 

2. DATA ACQUISITION AND INITIAL 
PROCESSING 

2.1. Observations and Image Reductions 
Ou r targets for opt i cal fo llow-up are UCHVCs selected 


Telesc ope to search for opti cal counterparts to a subset from Adams et al. (2013) as well as other ALFALFA 


of the Adams et al. (2013) UCHVCs and found a dis- sources that satisf y similar selection cr iteria but are not 


tant low surface brightness counterpart to an ALFALFA 


source (see also Sand et al. 2015 Adams et al. 2015b). 
This process - i.e., HI source ideiitihcation followed by 
optical imaging to identify a detectable stellar population 
and provide a distance estimate - has thus been shown 
to be a potentially important path for hnding more Local 
Volume dwarfs, for placing limits on how many nearby 
low-mass galaxies may remain undetected, and for un¬ 
derstanding how such previously-unknown dwarf galax¬ 
ies h t in with the rest of t he nearby galaxy population 
(e.g. McQuinn et al. |[2013 ). 


For the current project - follow-up imaging of 
UCHVCs with WIYN - we have devoted much of our 


AGC198606 
spite a radial 


Adams et al. 

|2015a and Table |2D 

ocitv below t 

Qe threshold lor incl 


in the Adams et al. (2013) sample (51 km s ^), the 


object has similar properties to Leo T’s HI component 
and is in fact near Leo T in both position and velocity 

^ The WIYN Observatory is a joint facility of the University of 
Wisconsin-Madison, Indiana University, the University of Missouri, 
and the National Optical Astronomy Observatory. 


part of the original Adams et al. (2013) catalog. We also 
add additional UCHVC candidates as the ALFALFA cat¬ 
alog expands beyond the initial 40% complete data re¬ 
lease. We observed the sources with the partially-filled 
One Degree Imager (pODI) on the WIYN 3.5-m telescope 
at Kitt Peak National Observatorj0 The pODI camera 
has a ^ 24' x 24' field of view and a pixel scale of 0.11" per 
pixel on WIYN. For each target, nine 300-second expo¬ 
sures were executed in a preset dither pattern in both the 
SDSS g and i filters, with some targets also observed in 
SDSS r. The majority of data collection for this project 
took place in March and April 2013, with additional ob¬ 
servations acquired in March, May, and October 2014. 
Future data collection will use an upgraded WIYN ODI 
camera with a 40' x 48' field of view. Observations and 
data analysis are ongoing; this paper presents the result 
of our efforts to establish the image processing meth¬ 
ods and analysis techniques for the project. The seeing, 

^ Kitt Peak National Observatory, part of the National Optical 
Astronomy Observatory, is operated by the Association of Uni¬ 
versities for Research in Astronomy (AURA) under a cooperative 
agreement with the National Science Foundation. 
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as characterized by the full width half maximum of the 
point spread function (FWHMPSF), in the images we 
have acquired so far ranges between 0.7" and 1.3" . The 
particular dataset presented in later sections is represen¬ 
tative of both the image quality and depth of the pODI 
observations we have obtained to date. 

Raw images were transferre d to the ODI Pipeline, Por¬ 


tal, and Ar chive (ODI-PPA Gopu et al. 2014 Young 
et al. |[20l4 1 at Indiana University, then processed using 


the OuickKeduce data reduction pipeline ( Kotulla |2014 1 
to remove the instrumental signature. The OuickKe- 
duce pipeline includes standard bias, dark, and flat 
corrections, as well as a pupil ghost correction, non¬ 
linearity corrections, cosmic-ray removal, and fringe re¬ 
moval where necessary. To account for image artifacts 
that remain after the QuickReduce processing, we con¬ 
structed and applied an illumination correction by cre¬ 
ating a dark-sky flat derived from the science observa¬ 
tions. To create the flat on a filter-by-fllter basis, we first 
masked all objects in the images, combined the masked 
images into a single image with a median algorithm, and 
then applied a smoothing kernel to the combined image 
to construct the final dark-sky flat field. We then divided 
each individual pODI image by the appropriate dark-sky 
flat field. This illumination correction was necessary in 
order to remove artifacts from mis-matched sky levels 
between dithered images. The final result is very flat im¬ 


ages over the full field-of-view of the camera (Janowiecki 


et al. 20151. The processed object images were then 


combined by scaling each image to a common flux level 
using Sloan Digita l Sky Survey (SDSS) DR9 catalog stars 
(Ahn et al. |2012 1 in the field, when available, to account 
for varying sky transparency. Additionally, we used the 
DR9 catalog stars to determine photometric solutions, 
and projected all individual images of the same field to a 
common pixel scale defined by the SDSS stars’ World Co¬ 
ordinate System (WCS) solution. Finally, the processed 
object images were stacked into a single combined image 
in each filter for each UCHVC held. 


2.2. Source Detection and Photometry 

A common set of steps is executed on each of the h- 
nal combined images. First, the images are cropped in 
order to remove the higher-noise edges of the dither pat¬ 
tern, leaving the central ^20' by 20' of the pointing. 
The mean standard deviation of the background level is 
measured in ~10—20 empty regions of the image. Using 
IRAF tasks daofind and phot, we identify all sources 
with peak counts at least 4.0 times the background noise 
level and measure their instrumental magnitudes. We 
explored the use of lower daofind thresholds but found 
that lower detection thresholds typically produced a large 
number of spurious sources without hnding any addi¬ 
tional genuine objects. We then mask out saturated stars 
(with diffraction spikes and/or bleed trails) and bright 
background galaxies and remove spurious sources with 
undefined magnitudes. To create a single list of sources 
appearing in both filters, we match sources in the g and 
i images and then measure the full width at half max¬ 
imum (FWHM) of the radial profile of each matched 
source. Obviously extended objects, with large FWHM 
values in a plot of FWHM vs. instrumental magnitude, 
are removed. The FWHM above which we remove ob¬ 
jects varies as a function of the instrumental magnitude. 


but a representative value is 1.5 times the mean stel¬ 
lar FWHM derived from bright stars in the image. An 
aperture correction is determined for each image and ap¬ 
plied to the measured instrumental magnitudes of the 
remaining sources. We convert these instrumental mag¬ 
nitudes to calibrated magnitudes by applying the zero 
points and color terms calculated from the SDSS stars 
in the images. Galactic reddening correctio ns are de¬ 
termined from Schlafly & Finkbeiner (201l|p and then 
applied to produce the hnal set of magnitudes and col¬ 
ors. At the end of this process, the pODI images with 
good seeing (FWHMPSF < 0.8”) typically yield >2000 
sources and a 5 (t limit on the brightness of a point source 
of ~25 magnitude in both g and i filters. 


3. DETECTION TECHNIQUE 
3.1. Color-magnitude filter 

Our objective is to search for the resolved stellar pop¬ 
ulation of the optical counterpart for each UCHVC we 
observe. We f ollow a si milar strategy to that outlined in 
Walsh et al. | (|2009|) and|Adams |(|2014l), making adjust¬ 
ments as appropriate for our data set. Ihe first step is to 
apply a filter in color-magnitude space to the calibrated 
photometry of the sources in each field. 

To constru ct a color-magnitude diagram (CMD) fil¬ 
ter, we select Girardi et al. (2004) isochrones for both 
old (8 < Cge/Cyr ^ 14) and young (8 < toge/M:yr 
< 14) stellar populations. Because isolated HI clouds 
are unlikely to be highly chemically enriched by the prod¬ 
ucts of stellar evolution, we use the metallicity range of 
Z = 0.0001 to Z = 0.0004 (1/200 to 1/50 of the solar 
value) subsets for each age group to further dehne the 
filter. The extent in color-magnitude space (i.e., in the 
Mi vs. g — i plane) covered by these isochrones is used 
to define the outer boundaries of the CMD filter. Figure 
[2 shows the isochrones selected to construct the hlter, as 
well as the final CMD filter we adopt. We note th at fo r 
some fields, such as the AGC198606 field (Section [4.2[ ), 
there is no evidence of bright, blue stars or ongoing star 
formation. In such cases we apply a CMD filter that in¬ 
cludes only the old stellar population, because there are 
simply no stars in the region of the CMD where young 
stars would appear. 

Before it is applied to the source photometry, the CMD 
filter is projected to a given distance within the Local 
Volume by adding the appropriate distance modulus to 
the Mi value of the original filter. We shift the filter 
so as to sample the range of distances over which we 
expect to be able to detect stellar counterparts to the 
UCHVCs: from 250 kpc (the closest distance at which 
we expect to find dark matter ha l os with a significant H I 
counterpart; |Adams et al. | (|2013|) , |Spekkens et al.| (|2014|) ) 
to 2.5 Mpc (the distance to the farthest observable Red 
Giant Branch star in a typical image from our data set). 
In terms of the distance modulus, we sample a range of 
m — M values from 22.0 to 27.0 in steps of 0.4 magnitude. 
When we find a potentially significant detection in the 
filtering/smoothing process, we re-sample the range of 
m — M values around which the detection occurs with 


^ Galactic extinction values were obtained from the NASA/IPAC 
Infrared Science Archive at http://irsa.ipac.calt6ch.edu/ 
applications/DUST/ 
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a more finely-spaced grid, in order to determine a more 
precise distance for the putative counterpart. 

The calibrated magnitude and color of each source is 
tested to determine whether the source falls within the 
CMD filter that corresponds to a specific distance. The 
source passes the test if its position in the i versus g — i 
plane falls within the boundaries of the CMD filter or 
if the la error bars (in either the magnitude or color 
direction) for the source overlap the filter boundaries. 
This has the effect of widening the CMD filter at fainter 
magnitudes. 


3.2. Spatial smoothing 

Once stars have been selected with the CMD filter at a 
given distance, we smooth the stars’ spatial distribution 
on the expected scale of a dwarf galaxy. This smoothing 
is achieved by first binning the remaining filtered stars 
into a grid with bin size ^ 7” x ~ 7”. We will refer to the 
bins in this grid as pixels, since in essence they represent 
a low-resolution image of the field. We then convolve 
each pixel in the grid with a Gaussian kernel with a size 
roughly equivalent to the spatial scale we expect for a 
nearby dwarf galaxy (FWHM ~ 2 arcmin). This spatial 
scale is comparable to the optical extent of the dwarf 
galaxies Le o T (half-light radius = 1.4 arcmin ; [Irwin 
et al. |2007l) and Leo P (optical extent ^ 90 arcsec;|Rhode 
et al. 2013). Pixels in the convol ved image (which ^ 


call A, following the convention of Walsh et al. 
with larger signal correspond to areas of the W1Y 


( 2 ^) 
N image 


with higher surface densities of stars that have passed the 
CMD filter criteria. We test whether these pixels contain 
genuine overdensities by calculating the mean signal level 
and its standard deviation over the entire A image. We 
then calculate a new image, S, such that 


S{x,y) 


A{x,y) - A 

-^(7 


( 1 ) 


where A and A^- are, respectively, the mean and standard 
deviation of the pixel values in A. This new image (S) 
effectively gives us the strength of the signal in A by mea¬ 
suring the number of standard deviations a given pixel 
is above or below the mean signal level, and provides an 
easy way to identify overdensities in the field. 

To evaluate the outcome of the filtering and smooth¬ 
ing process, we implement a technique that allows us 
to quantify the likelihood that an overdensity of a given 
strength (as defined in the S image) is a legitimate detec¬ 
tion of a possible optical counterpart rather than simply 
a stochastic variation. We use a Monte Carlo technique 
and generate 25000 samples of uniformly-distributed ran¬ 
dom points in a square 20' on a side, which is equal 
to the field-of-view of our cropped WIYN pODI images. 
The number of objects in each of the 25000 realizations 
is set equal to the number of point sources that pass the 
CMD filter criteria for a given field. We then perform 
the smoothing steps described above for each random 
realization and track the strength of the pixel with the 
highest signal level in the resultant S image. An example 
distribution of these peak values for a sample size of 400 
point sources is shown in Figure 

For each field and each application of the CMD filter, 
we can compare the peak overdensity value derived from 
the real, observed image data to the distribution of values 


generated from the random-sample data. In this way we 
can estimate the statistical significance of a given peak 
in the observed S image and assign a robust confidence 
level to detections. This is especially useful for images in 
which a marked overdensity is detected but which have 
no visually obvious optical counterpart. 


4. APPLICATION OF TECHNIQUE 
4.1. Leo P 

To illustrate our detection method, we apply it to im¬ 
ages of the galaxy Leo P, which was discovered by AL- 
FALEA, originally identified as a UCHVC based on its 


a dwarf galaxy in the Local Volume ( 

^ ^ -- r, - n 

Giovanelli et al. 

2013 Rhode et al. 12013 

AIcQuinn et al. 120131. We an- 

alyzed the original IRhoa 

e et al. 


images ot Leo R 


obtained with WIYN and Minimosaic in March 20IS_ 
Since these images were taken with broadband BVR fil- 
ters, the photometry measured by |Rhode et al. | (|2013|) 
was converted from BVR to q and i us i ng the stellar 


color transformations from Jester et al. (2005). After 
making slight adjustments to the smoothingand signifi¬ 
cance testing methods described in Section [31 to account 


for the difference in the field-of-view (Rhode et al. (20131 
analyzed a 9.6' by 4.8' portion of the Minimosaic point¬ 
ing) we applied our filtering and smoothing technique 
to the photometric measurements of Leo P. We utilized 
a CMD filter that included both young and old stellar 
populations as well as a filter that included only old stel¬ 
lar populations. Both filters produced highly significant 
detections; here we give the detailed results from apply¬ 
ing the old stellar population filter for better co mpa rison 
with the ACC 198606 results given in Section |4.2[ We 
searched over the range of distance modulus given above 
and found an extremely strong density peak, in all cases 
at least 8 -ct above the mean for the field, at each distance 
modulus. The peak occurred at the same spatial location 
in the field at all distances. In Figure [^ we show the re¬ 
sult of the filtering and smoothin g method applied at the 
measured distance of Leo P from McQuinn et al. (2013) 
(distance modulus 26.20; 1.738 Mpc), which is also the 
most significant detection of Leo P found with our search 
method. We begin by showing a map of detected sources 
in the image (the upper left panel of Figure [^, highlight¬ 
ing the sources selected by the CMD filtering process. We 
show the field CMD and filter in the upper right panel. 
The source map is then binned and smo othed according 
to the procedure described in Section [3(^ and is shown in 
the bottom left panel. Here we show a 2' diameter circle, 
centered on the location of the most significant density 
peak. In the bottom right panels, we show CMDs of the 
stars inside the circle (putative members of Leo P), and 
in a randomly placed reference circle. We see that the 
CMD of stars inside the circle clearly includes a vertical 
feature at {g — i)o = —1, associated with Leo P’s blue 
main sequence. This is an example of the kind of distin¬ 
guishing feature we might expect to see by examining the 
CMD of stars near the location of a stellar overdensity, 
if it is associated with a dwarf galaxy with active star 
formation. 


Images of Leo P obtained with pODI in April 2013 had poor 
seeing and could not be used for this analysis. 
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Figure 1. Construction of the co lor-magnitude filter we use to select stars consistent with an old or young stellar population. The colored 
points show|Girardi et al. | (|2004|| isochrones, and the shaded gray regions enclose our final filter boundaries. Dark gray is the old stellar 
population, and light gray is the young stellar population. Note that the old stellar population filter includes the region up to and including 
the tip of the Red Giant Branch, and does not include the Asymptotic Giant Branch region. 


Leo P is easily detected using the filtering and smooth¬ 
ing method. Figure|^shows the results of the significance 
testing corresponding to the peak at 1.738 Mpc and the 
CMD filter that includes only the old stellar population. 
The density peak associated with Leo P is an extremely 
strong detection at ~10.6-cr, with 123 stars passing the 
CMD filter, causing a density peak stronger than effec¬ 
tively 100% of the results from random distributions of 
stars spread over a field of the same size, with the same 
number of stars passing the CMD filter. These results 
verify that the methodology works robustly for the case 
of an object with a clear optical detection. We next ex¬ 
plore how the method performs on an object without an 
obvious optical counterpart. 


4.2. AGC198606 

AGC198606 is an ALFALFA-detected HI cloud with 
a dynamical mass of at least 3.5 x Kf'dMpc Mq, a ro¬ 
tational velocity of ~ 14 km s ~^, and a relatively l arge 
size on the sky of 23' x 16' (Adams et al. 2015a). A 
full list of Hl-derived properties can be found in 'Fable 
As explained in Section this s ource has a radial ve- 
locity too low to be include in the Adams et al. (2013) 
sample of UCHVCs, but its similarity to Leo 'L in terms 
of its HI characteristics, location, and velocity made it 
an interesting target for inclusion in our pODI imaging 
study. 

Observations of AGC198606 with WIYN and pODI 
took place in March 2013. Integration times of 2700 
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(T above local mean 


Figure 2. Distribution of peak a values for 25000 random samples 
with 400 points uniformly distributed across a 20^ by 20' field-of- 
view, which is well-fitted by a log-normal probability distribution 
function. 


seconds per filter in a nine-point dither pattern were ob¬ 
tained in both the g and i filters. Image reduction and 
photometric calibration were performed in the manner 
described in Section Errors on the photometric zero 
points and color terms ranged from 0.02 to 0.05 mag. 
Seeing was approximately 0.75" in both filters, with a 5- 
cr point source detection limit of 25.51 mag in g and 24.33 
mag in i. The final combined i-band image is shown in 
Figure along with 100" resolution HI contours from 


Adams et al. (2015a). 

We executed a series of artificial star tests to quantify 
the detection limits of the WIYN pODI images of the 
AGC198606 field. We began by quantifying the point 
spread function (PSF) in the g and i images using a set 
of bright, unsaturated stars. We added 300 artificial stars 
with magnitudes within 0.2 magnitude of a set value and 
the appropriate PSF to each of the g and i images. We 
then performed the same set of detection and photometry 
steps as performed on the original images and recorded 
the fraction of objects that were recovered in this process. 
We repeated these steps - i.e., adding 300 artificial stars, 
running the detection and photometry, and recording the 
results - in 0.2 magnitude intervals over a range of ~5 
magnitudes for each image. The 50% completeness levels 
are i = 24.1 mag and g = 25.1 mag. A dashed line in 
the upper right panel of Figure mark the i magnitude 
at which the convolved completeness (which takes into 
account the completeness level in both g and i) is 50% 
as a function oi g — i color. 

We applied the filtering and analysis steps described in 
earlier sections by shifting the CMD hlter for an old stel¬ 
lar population over the distance modulus range described 
in Section]^ from m—M = 22.0 mag to 27.0 mag in steps 
of 0.1 mag. We found three peaks of high significance in 
the density of point sources (p eaks that occur less than 
^ 10% of the time; see Section 3.2) that pass the CMD 
filter criteria near distance moduli of 22.4, 24.8, and 25.6, 
but these peaks were all located near the edge of the field, 
or outside the HI extent as shown in Figure Further 
inspection of these overdensities revealed that they were 
composed of faint extended sources not eliminated by the 
extended source cut and are probably background galaxy 
clusters. We also found an overdensity with ^ 90% con¬ 
fidence level near the center of the field at a distance 
modulus of 22.9. We more narrowly sampled the dis¬ 
tance modulus in steps of 0.01 mag around 22.9, and 
found that the peak with the highest significance is asso¬ 


ciated with a distance modulus of 22.89 (378 kpc) with 
total of 345 objects passing the CMD filter. A significant 
peak occurred at or near the same location with a simi¬ 
lar number of sources (344 - 360) with a confidence level 
of 67% or better between distance moduli of 22.86 (373 
kpc) and 22.97 (393 kpc). We also tested the effect of 
the size of the smoothing scale on the detected overden¬ 
sities. Larger smoothing scales do not reveal any unique 
strong overdensities in the image, while smaller smooth¬ 
ing scales typically find peak overdensities too small to 
be dwarf galaxies. 

Figure]^ shows the results of this analysis. The panels 
in this figure are identical to those in Figure except in 
the bottom left panel we indicate a 2' diameter circle cen¬ 
tered on the location of the most significant density peak 
in the image (magenta), the location of the HI centroid 
(black), and a reference circle with a randomly chosen 
location in the field (yellow). We select the stars inside 
the peak and reference circles and show CMDs for each 
subsample of stars in the bottom right panels of the fig¬ 
ure, allowing us to look for evidence of CMD features 
present only at the location of the density peak. We do 
not use the CMD for the entire field in this case since 
there are a much larger number of stars over a much 
larger area than in the Leo P field, and as such it would 
be difficult to distinguish features. The CMD centered 
on the density peak in this image exhibits a largely ver¬ 
tical morphology, parallel to the red giant branch from 
the isochrone-based hlter, while the CMD from the ref¬ 
erence circle does not include any objects that look like 
red giants at any distance. The dashed black line in the 
upper right panel indicates the 50% completeness level 
as a function of color and magnitude as described above. 
We also show typical instrumental magnitude and color 
errors as error bars along the right hand side of the up¬ 
per right hand panel. Typical errors were calculated by 
sorting stars into ten bins based on their apparent i mag¬ 
nitude and hnding the median instrumental magnitude 
and color error in each bin. Values for the instrumen¬ 
tal magnitude error are less than 0.01 for magnitudes 
brighter than i = 22 and increase to ^ 0.1 at a magni¬ 
tude of i ~ 24. Color errors reach the 0.01 level at a 
magnitude of i ~ 19.5 and the 0.1 level at i ~ 24. 

The density peak at m — M = 22.89 corresponding 
to 345 stars is ~ 4tT above the local mean for this field. 
Figure shows that the peak is higher than 92% of the 
peak values in a distributio n of 25000 random samples of 


the same size (see Section 


3.2|). Additionally, the peak 


locat ion is only ^2' from t he position of the HI centroid 
from Adams et al. (2015a|) (indicated by the black circle 
in Figure^ and still located within the HI disk. 

Figuref^shows a 5' by 5' cutout of the Aband image 
centered between the density peak (outlined with a ma¬ 
genta circle) and the HI centroid (marked with a black 
cross). There is no visually a pparent optical overd ensity 
in the image, as was noted in Adams et al. (2015a). The 
fact that some of the detected sources in the image are 
very faint and could be confused with unresolved back¬ 
ground galaxies, and the 92% significance level of the 
detection leave open the possibility that this is not a real 
detection of a stellar counterpart but is instead just the 
result of a chance superposition of sources with the right 
magnitudes and colors to pass through the CMD filter. 
However, the presence of some brighter stars along the 
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Figure 3. Results of the filtering and smoothing process for Leo P at a filter distance of 1.738 Mpc. Top left: sky positions of point 
sources relative to the corner of the field of view; small black dots are all sources detected in the image, red points are sources in the CMD 
filter. A 1 arcmin radius circle (magenta line) is centered on the detection peak. Another 1 arcmin radius circle (yellow line) is centered at 
a random point to provide a reference CMD. Top right: color-magnitude diagram for all point sources in the field; the CMD filter is the 
solid blue line, point colors as in the top left panel. Bottom left: the smoothed image S in units of standard deviations above or below 
the mean as described in Section [3^ where whiter pixels indicate higher density; the 1 arcmin radius circle (magenta) is centered on the 
highest signal pixel. Bottom right: color-magnitude diagrams for the stars inside the 1 arcmin radius magenta circle (left; colors as in other 
panels) and the stars in the yellow reference circle (right). 



Figure 4. Results of significance testing described in Section [3.2| 
on the detected overdensity in the Leo P field. Blue points are the 
histogram of density peak values for 25000 random distributions of 
123 points. The red line is the log-normal probability distribution 
function fit to the histogram. The dashed black l ine is the location 
of the Leo P detection at the |McQuinn et al. | (|2013|l distance of 
1.738 Mpc. - - 


red giant branch sequence, the location of the overden¬ 
sity near the center of the HI cloud, and the spatial and 
kinematic proximity of AGC198606 to Leo T, combined 
with the distance at which this overdensity has been de- 


tected (378 kpc, compared to 420 kpc for Leo T; Irwin et 


ar|[20^ , are compelling clues that this stellar overden- 


sity may indeed be a dwarf galaxy in the neighborhood 
of Leo T. 

We can use our data to place useful constraints on the 
total optical emission of the putative optical counterpart. 
We performed aperture photometry on the pODI images 
at the location of the detected stellar overdensity. Stars 
with i magnitudes brighter than 18 mag are brighter 
than the brightest star allowed in the CMD filter over 
the range of distance modulus used, and therefore are 
likely to be foreground Galactic stars and were masked 
in the images. We also masked red stars with {g — i)Z 2.0 
(guided by the CMD), and obvious background galax¬ 
ies in the vicinity of the overdensity and then measured 
the light within a 1.55' radius aperture centered on the 
location of the peak density. The size of the aperture 
was determined by reprojecting the 1.4' half-light radius 
of Leo T to the derived distance of the detected stellar 
overdensity at 378 kpc. 

After correcting for Galactic extinction and doubling 
the measured flux to account for only measuring in¬ 
side the half-light radius, we recover total magnitudes 
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Figure 5. WT YN nODT t-ha,nri imap- e of AGC198606, overlaid with 100” resolution HI column density contours derived from WSRT 
observations fseelAda ms et al. l2015all. The contour levels are fl, 2, 3, 4, SjxlO'^® atoms cm“^. The location of the HI centroid derived in 
|Adams et al. |(|20itia|| is marked witti a black cross. The field of view of this image is 20' by 20' and the orientation is N—up, E—left. The 
EWHMEbl of this image is approximately 0.75'' . We show the location of the detected overdensity in the field discussed in Section |4.2| 
with a 2' diameter magenta circle. 


of g = 19.23 ± 0.09 and i = 18.22 ± 0.09, with a, {g — i) 
color of 1.01. At a distance of 378 kpc as determined 
above, the i magnitude corresponds to an absolute mag¬ 
nitude Mi = —4.67. Note that these derived quantities 
contain light from unresolved sources in the AGC198606 
field, and so the luminosity and subsequent derived val¬ 
ues should be treated as upper limits. Using the rela¬ 
tionships between stellar mass-to-lig ht ratio and color s 

(| 200 ^, 


Bell et al. 


for u'g'r'i'z' passbands presented in 
we estimate the stellar M/L to be il.36, and at a dis- 
tance of 378 kpc, find an upper limit on the stellar mass 
of a potential optical counterpart to AGC19860 6 to be 


1.2 X 10'*^ M(t, . Gombined with the HI mass from Adams 
et al. (2015a), we find a Mhi/M steiiar = 42.5. 


We also derived lower limit values for the above quan¬ 
tities by adding the total flux in sources that lie within 
the GMD filter and are located inside the 1.55' radius 
aperture described above. This process gives us total 
magnitudes of * = 19.20 ± 0.09 and g = 20.18 ± 0.09, 
with a {g — i) color of 0.98, and an absolute magnitude 
Mi = —3.69. These values give a total stellar mass of 
4.6 X 10^ Mq and a Mm /M^tMa r = HO- The derived 
HI properties of AGGI98606 from Adams et al. ( 2015a| 
and our derived optical parameters are listed in Table [H 

5. DISCUSSION 

The putative optical detection of AGC 198606 is 
clearly provocative. For the purpose of the following dis- 
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Figure 6. Results of the filtering and smoothing process for AGC198606 at a filter distance of 378 kpc. Panels are the same as described 
as in Figure!^ In the upper right hand panel, typical magnitude and color errors are shown along the right hand side, and the dashed 
black line indicates the 50% completeness level. The 2' diameter magenta circle in the left hand panels corresponds to the location of the 
density peak discussed in the text. The 2' yellow circle is randomly placed and is used to select stars for a reference CMP, seen in the 
bottom right panel. Also plotted is a black circle indicating the position of the HI centroid reported in|Adams et al. |||2015aJ. The ~ 4 (t 
overdensity near the center of the field is a candidate for a stellar counterpart to the HI cloud. 


cussion we are going to assume that this optical detection 
is real, and that the stars identified in our images do 
in fact represent a sparse stellar population associated 
with the HI gas detected as a UCHVC by ALFALFA. 
We recognize that this interpretation is far from secure, 
but expect that future observations will be able to either 
confirm or refute the above hypothesis. 

The detection of such a small population of stars as¬ 
sociated with this HI cloud is consistent with models for 
what one might detect in a low-mass DM system that 
is at the boundary between systems with and without 
significant amounts of cold baryons. For exa mple, recent 
cosmologic al hydrodynamical simulations by |Ohorbe et| 
al. (2015) of isolated dark matter halos with masses of 
10^-5 - 10^° Mq result in galaxies with stellar masses of 
10"^ - 10® Mq and large Mg^g/M* ratios. The stellar mass 
estimates for AGC 198606 are at the l ower end of this 
range (10®'^ - 10^'^ Mq). Simulations by |Bovill fc Ricottl 


(2011) predict the existence of isolated ultra-taint dwarfs 


with similar properties. It is possible that AGC 198606 
represents exactly such a system as predicted by these 
simulations. 

To put the putative optical counterpart of AGC 198606 
into context, we compare its properties with several other 


dwarf galaxies in Table We include two gas-rich dwarfs 
(Leo T and Leo P) as well as three ultra-faint dwarfs 
(UFDs) that are gas-free companions to the Milky Way 
(Segue 1, Willman 1 and Segue 2) The char acteristics 


et al. 
from 
If 


for Leo P come from Rhode et al. (2013) and McQuinn 


2013), while those t or the other galaxies are taken 


cGonnachie 
indeed 


IS 


( 2012 ). 
real, 


the optical counterpart of 
AGC 198606 moves the detection of dwarf galaxies asso¬ 
ciated with UCHVCs into a new regime. The first discov¬ 
ery of an optical counterpa rt associated with a bona fide 
UCHVC wa s that of Leo P ( iGiovanelhet^aL TM^lRhodel 


et al. |2013| |Skillman et al.|2013| |Mc(4umn et al. 120131) 
but Leo P IS both substantially further away (1.73 Mpc) 
and larger (My = —9.4, Mhi = 9.3 x 10® M©) than 
AGC 198606. While the properties of Leo P are extreme 
(e.g., the lowest mass galaxy with current star formation, 
an abundance measurement consistent with the lowest 
metallicity systems known), AGC 198606 is even more 
extreme. The optical characteristics of AGC 198606 are 
more similar to those of th e UFDs discovered re cently in 
the halo of the Milky Way (|McConnachie |2012 and ref¬ 
erences therein). While possessing a more massive stellar 
population than the lowest mass UFDs listed in Table 
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Figure 7. WIYN pODI 2 -band image of AGC198606, shown in Figure [m cropped to a field of view of 5' by 5^ and centered near the 
density peak (marked with a 2' diameter magenta circle) described in the text. The orientation is N-up, E-left. The location of the HI 
centroid is marked with a small black cross. Stars that are included in the CMD filter are marked with red circles. 



Figure 8. Results of significance testing described in Section j3.2| 
on the detected overdensity in the AGC198606 field. Blue points 
are the histogram of density peak values for 25000 random distri¬ 
butions of 345 points. The red line is the log-normal probability 
distribution function fit to the histogram. The dashed black line is 
the location of the best AGC198606 detection at ~ 4cr. 


the total stellar mass in AGC 198606 is comparable to 
the U FDs recently reported in the Dark Energy Survey 
data (|The DES Collaboration et al.||2015|). 

When compared ^o lEs near neighbor Leo T, 
AGC 198606 exhibits the rather discrepant character¬ 
istics of having somewhere between ^1% and ~6% of 
the optical luminosity of Leo T but nearly double the 
HI mass. Not surprisingly, AGC 198606 has an extreme 
value for the HI gas to stellar mass ratio (between 45 
and 110). The HI distributions in the galaxies are rather 
dissimilar as well. For Leo T, the HI radius is 300 pc and 
the pe ak column density is meas ured to be 7 x 10^° atoms 
cm“^ (Ryan-Weber et al. 2008), while for AGC 198606 
the HI radius at the same column density level (2 x 10^® 
atoms cm“^) is ^ 2 times as large (600 pc), the full HI 
extent is ~ 4 times larger (r = 1250 pc), and the peak 
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Table 1 

Properties of AGC198606 and Possible Optical Counterpart 


Property 

Value 

HI properties 

lo: 

09:30:02.5 

Dec 

-1-16:38:08 

Mhi 

3.5 X 106^2^ Mq 

dni 

11' ± 1' 

a X b 

to 

CO 

X 

rni 

3.3 dMpc 

^dyn 

3.5 X 10^ dMpc 

Nhi 

6 X 10^® atoms cm~^ 


Estimated Optical Properties 


bright limit faint limit 

Distance 

378 kpc 

a 

19.23 di 0.09 mag 20.18 di 0.10 mag 

i 

18.22 ± 0.09 mag 19.20 ± 0.17 mag 

Mi 

-4.67 ± 0.09 mag -3.69 ± 0.17 mag 

Stellar luminosity 

5.0 X IO^Lq 2.0 X 10^ Lq 

Stellar mass-to-light ratio 2.36 2.27 

Stellar mass 

1.2 X lO'^M© 4.6 X lO^M© 

Mni/MsteUar 

42.5 110 


_ N ote. — HI properties reproduced from Table 1 of [Adams et| 

jal. |||2015a|| . Stellar mass-to-l ight ratios have been estiniated from 
relations in [Bell et al. | l (2003[ l 


column density is a factor of ^ 10 lower (6 x 10^® atoms 
cm“^). We note that the beam sizes used for measuring 
the peak column densities for the two galaxies are similar 
(- 45" for Leo T and ~ 60" for ACC 198606) but the 
HI extent of AGC 198606 is measured using the low res¬ 
olution 210") data; these mismatches should be taken 
into account when considering the numbers above. 

Prior to the discovery of A GC 198606, Leo T ap peared 
to be fairly isolated in space ( McConnachie |2012[ ). With 
a projected angular separation ot 1.2"", the minimum 
physical separation of the pair is only 8.4 kpc (adopt¬ 
ing the average distance of the two galaxies of 400 kpc). 
Hence, Leo T and AGC 198606 may be separated by only 
a modest number of HI radii (~6 x ihi for AGC 198606). 
More than likely their current physical separation is 
larger, perhaps 30-40 kpc. However, given their simi¬ 
lar distances and small velocity difference, it would seem 
plausible, even likely, that the two galaxies have had one 
or more close encounters in the past. It is tempting to 
associate the fairly rece nt episode of star formation in 


Leo T (^ 200 Myr ago; Irwin et al. 2007 de Jong et 


al. 2008) with such an encounter. The apparent lack 


ot any obvious tidal HI bri dging between the two galax¬ 
ies in the ALFALFA data (Adams et al. 2015aI would 
suggest that any such encounter was not a strong one. 
The presence of foreground Galactic emission, however, 
limits the constraints that can be put on any connecting 
emission. 

It is even possible that Leo T or its “friend” has a 
tidal origin and originated from the other. Given the 
more substantial stellar population associated with Leo 
T, the more natural interpretation would be that such 
a tidal encounter with a third system might have pulled 
roughly two-thirds of the HI gas out of Leo T to form 
AGC 198606. If so, then it would appear that any tidal 
model would need to be finely tuned to explain the ob¬ 
served system. That is, it is hard to imagine a tidal event 
that could remove more than half of the gas from such a 


low-mass system without disrupting the system entirely. 
The absence of an obvious perturber is another problem 
with this scenario, as is the lack of tidal HI mentioned 
above. We do not favor a tidal origin for AGC 198606, 
although such a scenario cannot be ruled out. 

As mentioned above, the optical counterpart to 
AGC 198606 is much more similar to the UFDs than to 
the more gas-rich systems like Leo T or Leo P. It has been 
assumed that some UFDs represent the remnant stel¬ 
lar populations of heavily s tripped dwarf galaxie s that 


were accreted by the MW (Willman et al. 


2006 Mar¬ 


tin et al. 2007). The existence ot the extremely sparse 


stellar population in AGC 198606 raises the possibility 
that the origin of these UFDs may have been with sub¬ 
stantially less massive systems than has been previously 
supposed. That is, UFDs may have had progenitors that 
looked more like AGC 198606 than more massive dwarf 
spheroidals. In this case, far less stripping would be re¬ 
quired to account for the low stellar masses of some cur¬ 
rent UFDs. While there is no reason to doubt the va¬ 
lidity of the stripping scenario for the creation of some 
of the UFDs, the existence of systems like AGC 198606 
provides another creation path for this class of objects. 

It seems clear that additional observational work is 
needed before any optical detection of a stellar com¬ 
ponent to AGC 198606 can be considered secure. The 
sparse nature of the putative stellar population will make 
a definitive detection challenging. Two logical steps for¬ 
ward would be (1) deeper ground-based imaging and (2) 
spectroscopy of the brighter RGB stars detected in our 
current data. Our current pODI images are already fairly 
deep and are of good image quality. At the location of 
the RGB in our GMD the pODI data are 50% complete 
at * = 24.1. However, the location of the red clump at 
i ^ 23.5 falls where both field contamination and in¬ 
creasing photometric errors combine to make a defini¬ 
tive detection of the red clump stars uncertain. Deeper 
data with reduced photometric errors would help to im¬ 
prove this situation. Radial velocity measurements of the 
handful of possible RGB stars in AGC 198606 have the 
potential of confirming or rejecting their possible asso¬ 
ciation with the HI gas. Such observation s were useful 
in the case of Leo T (Simon & Geha 2007) and are well 
within the reach of current ground-based spectroscopic 
capabilities. 

Our group has obtained optical imaging data of sim¬ 
ilar degthUo^jn^y additional UCHVCs selected from 
(|2013|) catalog and from subsequent 
WEFA data by members of our team. 


the Adams et al. 
searches ot the AL 
Initial visual searches of our current imaging data reveal 
no other optical counterparts as obvious as Leo P. We 
are currently analyzing the data set using the automated 
method described here, with the expectation of either de¬ 
tecting additional sources similar to AGG 198606 or set¬ 
ting upper limits on the mass of any stellar component 
that might be present. 


6. SUMMARY 

We have presented a processing and filtering technique 
designed to detect possible optical counterparts of ultra¬ 
compact high velocity HI clouds identified by the AL¬ 
FALFA survey. Our filtering technique relies on con¬ 
structing a region in color-magnitude space consistent 
with the theoretical and observational properties of other 
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Table 2 

Properties of Dwarf Galaxies 


Galaxy 

Distance 

My 

re 

Mhi 

M* 

Mhi/M i, 

“Friend of Leo T” (AGC198606) 
Leo T 

Leo P 

Willman 1 

Segue 1 

Segue 2 

378 kpc 

409 to 420 kpc 
1720 kpc 

38 kpc 

23 kpc 

35 kpc 

-4.7 to -3.7 
-8.0 
-9.4 
-2.7 
-1.5 
-2.5 

120 pc 

25 pc 

29 pc 

35 pc 

5.0 X 10® 

2.8 X 10® 

9.3 X 10® 

4.6 - 12.0 X 10® 

~ 2 X 10® 

5.7 X 10® 

42 to 110 

~ 1.4 
~ 1.6 


Local Group compact dwarf galaxies, which we then ap¬ 
ply to the photometry of point sources in WIYN pODI 
images in order to detect overdensities in the point source 
maps. This method provides an effective way to de¬ 
tect stellar overdensities given appropriately deep im¬ 
ages. Applying the technique to photometric observa¬ 
tions of Leo P, a recently discovered star-forming dwarf 
galaxy, we recover an extremely significant overdensity 
in the image, implying that less obvious overdensities 
should also be detectable via our filtering and smoothing 
process. 

We find compelling evidence for the presence of a com¬ 
pact dwarf galaxy near the location of AGC198606, a 
UCHVC in the vicinity of Leo T, the faintest dwarf 
galaxy in the Local Group with recent star formation 
and HI gas. Our inferred properties for a potential dwarf 
galaxy give it a diameter of ~ 2 arcmin at a distance of 
^ 378 kpc, only 42 kpc less distant than Leo T. An inter¬ 
esting candidate for follow-up study, the possible optical 
counterpart to the UCHVC AGC198606 is a good ex¬ 
ample of the capabilities of our processing and filtering 
methods. We are just beginning to apply the methods 
to the WIYN pODI imaging data we have in hand and 
hope to report on additional detections, or at least to set 
limits on the stellar populations that may be associated 
with the ALFALFA UCHVC sources we have targeted. 

We are grateful to the staff at the WIYN Observatory, 
Kitt Peak, and NOAO for their help during our various 
WIYN pODI observing runs. We also thank the ODI- 
PPA support scientists and programmers at WIYN, lU, 
and NOAO for their support and assistance with the re¬ 
duction and analysis of the WIYN pODI data. The au¬ 
thors acknowledge the work of the entire ALFALFA col¬ 
laboration for their efforts in support of the survey. We 
thank Tom Oosterloo for useful discussions, especially 
about the HI observations of this object. We would also 
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AST-1107390 to R.G. and M.P.H. and by grants from 
the Brinson Foundation. J.M.C. is supported by NSF 
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